Multiphoton above threshold effects in strong-field fragmentation 
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We present a study of multiphoton dissociative ionization from molecules. By solving the time- 
dependent Schrodinger equation for and projecting the solution onto double continuum scattering 
states, we observe the correlated electron-nuclear ionization dynamics in detail. We show — for the 
first time — how multiphoton structure prevails as long as the energies of all fragments are accounted 
for. Our current work provides a new avenue to analyze strong-field fragmentation that leads to a 
deeper understanding of the correlated molecular dynamics. 
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Despite more than 20 years of scrutiny, strong-field dis- 
sociative ionization of molecules is still not completely 
understood. Understanding this process, however, would 
provide insight into how the energy deposited in the 
molecule by an intense laser pulse is shared between the 
electrons and the nuclei via their correlated motion. A 
large part of the challenge in investigating this process 
is that the dynamics of an ionized electron is not eas- 
ily treated by the usual Born-Oppenheimer (BO) ap- 
proximation. Complicating the issue further is the fact 
that the final state lies in at least a double continuum, 
likely Coulombic, comprised of the free nuclear and elec- 
tronic motion which raises fundamental questions about 
the analysis. 

Because ab initio calculations of dissociative ionization 
require going beyond the BO approximation, the vast 
majority of intense field dissociative ionization calcula- 
tions have been carried out for the simplest molecule: 
(see, for example, Refs. [ij-Q)- Even for this system, how- 
ever, full-dimensional solutions of the time-dependent 
Schrodinger equation (TDSE) have not yet been obtained 
in an intense field as they pose an immense computa- 
tional challenge. Consequently, the dissociative ioniza- 
tion calculations that have been done simplified the prob- 
lem even further through ad hoc reductions of the dimen- 
sionality or other severe approximations. 

Nevertheless, intense field processes identified via Hj' 
studies — even reduced dimensionality ones — are now 
understood to occur in other systems as well. Dissocia- 
tion via bond softening and above threshold dissociation 
(ATD), for instance, were first identified in and have 
now been seen in other systems like 00^+ [1], Oj [3, 
NaJ [6], and 0. 

Ionization mechanisms such as charge-resonance en- 
hanced ionization (CREI) have also been identified in 
@, 01 and applied to other systems [13, HI- More 
recently, two new models have been proposed to explain 
unexpected structure measured in the nuclear kinetic en- 
ergy release (KER) spectrum following ionization of H J . 
In one model, the interference of two dissociation path- 



ways leads to the modulation of the KER spectrum via 

CREI [mil. The other model, named above thresh- 
old Coulomb explosion ( ATCE) , is based on the Floquet 
potentials obtained by dressing not only the field-free 
BO potentials but also the 1/R Coulomb explosion curve 
with photons from the laser field (see Fig. 2]) [l3|. The 
observed structure in the KER spectrum is thus due to 
absorption of different numbers of photons [3, HI ■ The 
predictions of the models deviate for low intensities and 
there is no consensus on which is correct. So, in addi- 
tion to answering fundamental questions about electron- 
nuclear correlations, TDSE solutions — along with an 
accurate way to analyze them — are needed to resolve 
this controversy. 

Key to understanding electron-nuclear correlations is 
the identification of an appropriate physical observable. 
It has long been recognized [l6| that a particularly useful 
observable for this purpose is the joint energy spectrum 
(JES). For instance, experiments at the Advanced Light 
Source have shown how the simultaneous measurement of 
the KER and the energy of a freed electron can be used to 
extract details about the single-photon- induced breakup 
of molecules such as C2H2 [l3| and CO [S]- Joint en- 
ergy spectra are also useful beyond electron-nuclear dy- 
namics. In fact, they provide insight whenever the fi- 
nal state involves three or more fragments. An applica- 
tion to the double ionization of He in an intense laser 
field showed, for example, that the joint electron en- 
ergy spectrum would reveal a clear separation of sequen- 
tial from non-sequential processes 1 1911 and could thus 
help answer an outstanding question |20| in intense field 
physics. While not exactly a JES, the goal of the well- 
known Dalitz plots [il'l is very similar. Yet, even with so 
much evidence of their utility, JES following intense field 
dissociative ionization have not been studied. 

Our goal in this Letter is to show that the JES does 
indeed provide considerable insight into strong-field dis- 
sociative ionization. First and foremost, it shows clear 
multiphoton structure reminiscent of — but distinct from 
— above threshold ionization (ATI) ATD and 
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ATCE that was not previously seen or anticipated. But, 
because the JES clearly shows how the energy is shared, 
it can also help answer the outstanding question of how 
much energy is taken by the electron during strong-field 
ionization. Finally, the JES may well provide an alter- 
native and detailed avenue to time-dependent imaging of 
molecular processes. 

While our conclusions generalize to any molecule, we 
will illustrate them here for dissociative ionization of H J . 
The simplicity of this system removes the theoretical ap- 
proximations that would otherwise be present for multi- 
electron or polyatomic species, providing an unambigu- 
ous demonstration of our ideas. For this same reason, 
we will also consider a reduced-dimension model of 
to allow for an essentially exact numerical solution of the 
TDSE. Since our results are based on energy distribu- 
tions and depend primarily only on energy conservation, 
they are not sensitive to the dimensionality, making the 
reduced-dimension the most transparent example for 
our purposes. 

Specifically, we consider a linearly polarized laser field 
with the nuclei aligned along the polarization axis, in- 
cluding just the internuclear separation R and the elec- 
tronic coordinate x in the direction of the laser polar- 
ization, measured with respect to the nuclear center-of- 
mass. The length gauge Schrodinger equation then reads 
(atomic units are used throughout), 

i^^^{R,x,t) = [Ht,+H, + £{t)x]^iR,x,t), (1) 
ot 



where 
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with XA,B ^ X ± R/2. We have softened the Coulomb 
singularity with the function a(i?), chosen to reproduce 
the full-dimensional lag BO potential — and thus the 
ionization potential [24,-i26i] . We use a laser electric field 
of the form 



£{t) ^ £q siu^iTTt/T) cos(wt), < t < T 



(4) 



with T — N{2t:/lu) in which N is the number of cycles 
(tfwhm — 0.364r). Results in this paper use N = 10. 

We solve Eq. ([1]) using the finite element discrete vari- 
able representation [27H30j . The {R, x) plane is divided 
into inner, |a;| < 50 a.u. and i? < 15 a.u., and outer, 
50 a.u. < |a;| < 1500 a.u. and 15 a.u. < i? < 50 a.u., re- 
gions. Within each region, R and x each have a uniform 
element distribution with the element size based on the 
minimum expected de Broglie wavelength. The abrupt 
change in element size at the boundary between regions 
is eliminated by repeated application of three-point av- 
eraging. We propagate the solution using the short-time 



evolution operator, 

X,t + S)^ ^-t[H^+H^+£it+5/2)^]S^^j^^ 

evaluated with the Lanczos algorithm [3l| . This com- 
bination of techniques produces very accurate solutions 
across many orders of magnitude and, based on testing, 
gives an accuracy of at least two significant digits. 

A critical aspect of calculating the electron- nuclear en- 
ergy spectrum is the careful separation of the p + p + e~ 
double continuum from the p + H single continuum. We 
accomplish this separation by calculating the double con- 
tinuum wavefunction -0^^ ^ within the BO approxima- 
tion. This approach guarantees their orthogonality to the 
bound electronic states, and it incorporates the Coulomb 
interactions. The probability of observing a KER of i?N 
and an electron energy of E^. is thus 



dE^dEe 



\^{t))\ 



(6) 



where g, u denotes gerade and ungerade symmetry, re- 
spectively. Explicitly, our approximation to the double 
continuum scattering states is 



(7) 



with X-En s-i^d (jf^ energy normalized scattering states 
satisfying 



H^XE^iR) — -Enx_Bn(-R), 
H^4>%{R-x) = E,(j,%{R-x) 



(8) 
(9) 



which are each solved using a variational R-matrix for- 
mulation [3^ . 

Since this analysis method has not been applied be- 
fore to strong-field processes, we emphasize its attrac- 
tive features beyond the rigorous separation of the dou- 
ble continuum it provides. Because it is based on using 
the energy eigenstates, the spectrum can be calculated 
as soon as the pulse is negligible. So, while the result- 
ing spectra are consistent with existing approaches such 
as a simple Fourier transform 33| and the scaled coor- 



dinate approach 3J|, our method is much cheaper com- 
putationally as we need not propagate the solution to 
macroscopic times — the latter being necessary both to 
make the Fourier transform a good approximation for a 
Coulombic system 33| and to improve the approximate 
identification of the double continuum based on spatial 
position. 

We solved Eq. ^ for wavelengths from 230 to 650 nm, 
for intensities from 10^'^ to lO^'' W/cm^, and for a number 
of initial vibrational states. As measured by the Keldysh 
parameter [35| . all cases are weakly within the multi- 
photon regime: neglecting nuclear motion and using an 
ionization potential of 1.1 a.u., the Keldysh parameter 
lies in the range 2-8. 
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FIG. 1: (Color online) JES from Eq. (O for exposed to a 
400 nm, 8.8x10^^ W/cm^ laser pulse starting from (a) u = 0, 
(b) V = 2, (c) V = 7, and (d) w = 9. The thick black curves 
in the top and side panels show the ATI and KER spectra, 
respectively, from integration of the JES. For comparison, 
frozen nuclei and reflection results (see text) are shown with 
thin grey lines. 




FIG. 2; (Color online) The JES from Eq. Q starting from 
w = 9 at 5.6 X 10^^ W/cm^ as a function of wavelength: (a) 
650 nm, (b) 506 nm, (c) 400 nm, and (d) 363 nm. 



Figure [T] shows JES for exposed to a 400 nm laser 
pulse, starting from v—0, 2, 7, and 9 of HJ(1(7). 

A characteristic feature of these density plots is the 
maxima along lines given by Ey + nhuj = E-t^ + E^ + Up. 
Here, Up is the ponderomotive energy of the electron due 
to the presence of the laser field |47|; and E^, the en- 
ergy of the initial state relative to the p + p + e~ thresh- 
old. In other vifords, these plots are evidence of multipho- 
ton absorption in the molecular double continuum. This 
conclusion has been confirmed via calculations at other 
wavelengths (see Fig. [2]). 

Interestingly, these JES reveal much more than multi- 
photon features. For instance, in the present case they re- 
flect strong electron-nuclear correlation. To understand 
this, we first note that — as with any homonuclear di- 
atomic molecule — the laser field couples only to the 



relative motion of the electron and the nuclear center of 
mass. Thus, energy can be transferred into the relative 
nuclear motion only through interactions with the elec- 
tron. Yet, for a given number of photons absorbed, there 
is a nonzero probability that the relative motion of the 
nuclei takes most of the energy, leaving the electron with 
almost no asymptotic kinetic energy. This behavior is re- 
vealed in Fig.[T]by the broad distribution in i^N at Ee ^ 
which indicates efficient energy transfer from the electron 
to the nuclei. This strong correlation could be important 
in other molecules, thus raising questions about neglect- 
ing nuclear motion even for short laser pulses. 

In order to understand the additional structure ap- 
parent in Figs. [T] and [2l we consider the strong-field ap- 
proximation (SFA) [36| for the case of dissociative ion- 
ization [ST*! and approximate the final state of the nuclei 
with a plane wave. We note, in passing, that a somewhat 
similar extension of the SFA has been made to study the 
correlated dynamics of multi-electron atoms [SS*] . For the 
present case, one can show that within the BO approxi- 
mation, the SFA gives [47] 



d'Ps 



FA 



dE^dE^ 



oc |x«(v^TOpi?N) 



(10) 



where Xy and 0o are the Fourier transforms of the initial 
vibrational and electronic wave functions, respectively. 
The presence of \xv\^ explains the origin of the features 
at constant E-^^. The factor |i^oP similarly shows that fea- 
tures at constant Ef. are also possible. Together, these 
factors show that the JES can be used to image the to- 
tal electron-nuclear initial state. The factor F^^ ^ in 
Eq. ([T0|) contains all of the laser parameter dependence 
and thus accounts for the multiphoton lines in the JES. 

In all previous studies of this system (see, for example, 
Refs. [i,ll2, 15, l^li^l) — and multiphoton dissociative 
ionization more generally — it was the nuclear KER spec- 
trum that was presented and not the JES. Moreover, a 
substantial fraction of experimental and theoretical work 
still relies on the simple refiection method either to pro- 
duce the KER spectrum from a nuclear wave packet upon 
ionization or to deduce the nuclear i?-distribution from a 
measured KER spectrum 41-41]. Since the pulse lengths 



we consider are all short compared to the vibrational 
periods, the reflection method amounts to the mapping 
dP/dEN K \xvil/EN)\yEl,. 

With our accurate analysis of the double continuum, 
we can test this simple method for our model HJ. Our 
KER spectra, obtained by integrating the joint spectra 
over Ee , are thus shown in the side panels of Fig. [1] along 
with the reflection method spectra. Although the reflec- 
tion method accounts for some of the gross features of the 
KER, key features are not reproduced. This failure sug- 
gests the violation of at least one of its foundational as- 
sumptions: (i) that the ionized electron is fast and leaves 
behind two bare protons that subsequently Coulomb ex- 
plode, (ii) that the nuclei have zero kinetic energy to 
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start this explosion, and (iii) that the ionization rate is 
independent of R. 

Just as the KER spectrum has been used to study nu- 
clear dynamics, much has been learned about electronic 
dynamics by looking at the electron's ATI spectrum. The 
existence of the ATI peaks, for instance, reflects the peri- 
odic launch of electron wave packets at each half-cycle of 
the laser pulse. Additional features in the spectrum arise 
from electronic structure and intra-cycle dynamics (45| . 
The top panels of Figs. [1] and [2] show the results of in- 
tegrating the JES over i?N, and, interestingly, the ATI 
peaks survive. Their survival is a consequence of the fact 
that the JES peak sharply at constant i?N- 

For comparison, we have calculated the ATI spectra 
by solving the TDSE for our model with the nu- 
clei at a fixed distance corresponding to the maximum 
of \xv{R)\^- The results are shown in the top panels of 
Fig. [1] as well. We note that these ATI peaks are both 
more pronounced and shifted compared to those obtained 
from the full calculation, underscoring the importance of 
nuclear motion even for very short laser pulses. 

From Figs. □ and H it is clear that the KER and ATI 
spectra — being projections — obscure much of the in- 
formation visible in the underlying JES. We thus suggest 
two different and complementary spectra for studying the 
dynamics: the total energy and energy sharing spectra. 

Since the multiphoton peaks occur along energy con- 
servation lines, it is natural to integrate along lines of con- 
stant energy to obtain the total energy spectrum. Car- 
rying out the integration for the JES of Fig. [TJ we obtain 
the results shown in Fig. [3ja). The opening of lower-n 
photon channels with increasing v apparent in the fig- 
ure is readily predicted from the ionization-supplemented 
Floquet curves in Fig. In particular, the lowest peak 
visible for a given v is determined from Fig. 3] by the 
lowest-n ionization curve that crosses the Ig — Ow curve 
at a distance where \xv{R)\'^ is non- negligible. For in- 
stance, Fig. m shows that n=6 crosses Ig — Ouj, but it 
does so where \xv{R)\'^ is negligible for all v in Fig. |31 
Consequently, the lowest peaks visible have n=7. But, 
for v=0 and 2, even the n=7 crossing is not accessible, so 
their lowest n is 9 and 8, respectively. This finding lends 
support the validity of the ATCE model which is based 
on such a figure [ij, [3 ■ 

It is important to emphasize that we expect a figure 
like Fig.|3Ua) for any system so long as the energies for all 
the fragmentation degrees of freedom are included in the 
total energy. This expectation follows from the fact that 
the occurrence of multiphoton features depends essen- 
tially only on energy conservation. By extension, we thus 
expect to see multiphoton features like those in Figs. [1] 
and m in the JES for any system when all the fragmen- 
tation degrees of freedom are accounted for. Moreover, 
the same conclusion is reached by generalizing the SFA 
result in Eq. (|10p to other systems. 

These multiphoton features thus provide a natural bin- 
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FIG. 3: (Color online) (a) The total energy {E = En + Ee) 
spectra for the JES in Fig. [T] (b) Cuts of the JES from 
Fig[ljc) demonstrate the energy sharing between the electron 
and nuclei starting from v = 7 for 7- {E — 0.3 a.u.), 8- {E = 
0.4 a.u.), and 9-photon {E — 0.5 a.u.) absorption. 
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FIG. 4: (Color online) Dressed bound and ionization thresh- 
old potentials at 400 nm for our ID Hj" model. The horizontal 
lines indicate the energies of the states v = 0,2,7, and 9. 



ning, suggesting that we study the energy sharing spec- 
trum for each photon peak. Figure [3l^b) shows precisely 
this for the v—7 state. Note that the structure from the 
vibrational state is still clearly visible. Since the spec- 
trum for each n peaks at the same E-^, the nuclei tend 
to take the same amount of energy — the minimum pre- 
dicted from Fig. 2] — while the electron tends to take all 
energy from the excess photons, explaining the survival 
of the peaks in the ATI spectrum. 

The tendency we note for the electron to take most 
of the excess photon energy may generalize to other sys- 
tems, but further investigation is needed. What is clear, 
is that a straightforward extension of the SFA result in 
Eq. (ITOl) will apply to other systems implying that the 
Fourier transform of the initial state will be imprinted 
on each energy sharing peak. 

In this Letter, we have highlighted just some of the in- 
sights possible from the JES. By illustrating them with 
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a physically transparent reduced-dimension model, 
we have obtained some initial answers to several ques- 
tions for dissociative ionization such as how the energy is 
shared between the electron and nuclei and how the JES 
might be used for imaging the total — electron plus nu- 
clear — initial state. Further study is certainly in order, 
and this Letter suggests many possible directions. For in- 
stance, while our SFA extension already shows the imag- 
ing possibilities, extending it a little further hints that the 
JES might also be a sensitive indicator of the breakdown 
of the BO approximation in the initial state. In addition, 
we believe that studying the JES for IR-pump and XUV- 
probe schemes will provide the additional data needed to 
establish the validity of the ATCE. As it happens, strong- 
field experimentalists are currently tantalizingly close to 
being able to measure three particles in coincidence with 
sufhcient statistics to produce a JES 44, 4^. 
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I. STRONG-FIELD APPROXIMATION 



In this supplementary material, we discuss the extension of the SFA to the case of disso- 
ciative ionization. Our approach closely follows the methodology of [37], but we approximate 
the final state of the nuclei with a plane wave [see Eq. (jl]) below] instead of a Coulomb wave. 
We do this to consistently treat the electron and nuclei at the same level, since our TDSE 
calculations indicate that the effects of the Coulomb potential on the nuclei and the electron 
are equally important. More importantly, doing so allows us to derive the simple formula 
given in Eq. (10) of the main text. This simple formula for the JES is consistent with 
the TDSE result. Further, it highlights the imaging application of the JES and potentially 
provides a much less computationally demanding means of computing the JES. 

Within the S-matrix formalism, the joint energy distribution can be written as 

d^PsFA 



^2n . r 2 



oc 



dE^dE^ /-f \k^\\ke 



(1) 



where the sum is over four terms, corresponding to the positive and negative values of the 
final momenta /cn and k^. leading to the same energy. In each term, VExt represents the laser 
molecule interaction, i.e., E{t)x in the length gauge and A{t)d/dx + A'^{t)/2 in the velocity 
gauge, with A{t) = — /* dt'£{t') the vector potential. 
As the initial state, we use 

ij,{R,x,t)=Xv{R)Mx)e-'^^\ (2) 

which presupposes the validity of the Born-Oppenheimer (BO) approximation and that the 
electronic wave function may be approximated near the nuclear equilibrium position as an 
i?-independent function (note that improving this assumption by using linear combination 

— 1/2 

of atomic orbitals will lead to the same result within the order of nip ). These assumptions 
are justified in this work by the fact that all considered laser pulses have a duration short 
compared to the vibrational period. Consistent with the SFA, we employ a product of a 
plane wave (PW) and a Volkov wave for the final state 

7A/(fcN, fee; R, X, t) = xl^iR, t)C°'(^. t). (3) 



2 



We proceed using the velocity gauge, where 

,PW/ 



ik^x-i^ f{k^+A{t')fdt' 



Inserting these into Eq. ([T]), we have 



27r 

1 ^ik^R-i-±-klt 



(4) 
(5) 



SFA 



oc 



^i[{E^+E^-E^)t+f[A\t')/2+keA{t')\dt'] 



dt 



(6) 



with i?e = k1/2, = k'^/rrip, and Xv,4>o denoting the Fourier transform of the respective 
wave functions. 

We note that in the length gauge, the Fourier transform of the electronic state is evaluated 
at momenta shifted by the vector potential at each time and thus does not separate outside 
the integral. Additionally, the velocity gauge interaction is replaced by that of length gauge. 
However, by evaluating the time integral using the saddle-point approximation, one will 
arrive at the same qualitative statements as we derive below using the velocity gauge. 

Our current numerical investigations are in the regime where multiphoton structure is 
displayed. We may thus assume that the vector potential has the simple form 



A{t) = Aocos{ut). 



(7) 



We may then derive an explicit expression for the joint energy distribution: 



dE^dEe \k^\\ke\ 

^ u ' iujj J 

^ |x.(v^^^V^)|'|0o(\/2^)l' 



A{t)ke + 



dt 



yj2mpE^Ee 

^(_1)"J„ /-/2^^^_A|\ I ^\^^E,+E^+Al/i-E^- 



-nuj)t\ 



A{t)y/2Ee + 



A\t) 



dt 
(8) 
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where 
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n 




2 



dt 



2 



(9) 



with Jn a generahzed Bessel function [36]. To get rid of the summation over nuclear and 
electron momenta in Eq. ([8]), we exploited the symmetry of the initial vibrational and elec- 
tronic wave functions along with the following properties of the generalized Bessel function 
and the vector potential: J„(— = (— l)"J„(M,f) and Ait + n/oj) = —A{t). We also 
exploited the relation Jn{u, —v) = (— 1)" J_n(M, f ) to change the sign of n in the phase of 
the time integral. 

In Eq. ([9]), the time integral picks out the energy conservation lines as defined by 



where n corresponds to n-photon absorption and Up = A^/A represents the ponderomotive 
energy shift that the electron acquires due to the oscillating field. The generalized Bessel 
function accounts for the observed laser parameter dependence of the individual multiphoton 
peaks. We point out that for larger molecules, the Fourier transform in Eq. ([8]) is replaced 
by the Fourier transform of the total wave function and that the sum E-^ + in the phase 
of the time integral is replaced by the total energy. Additionally, the first argument of the 
generalized Bessel function and the laser- molecule interaction part, VExt, of the time integral 
will have additional terms for any particle that couples directly to the laser field. From these 
observations, it immediately follows that our findings generalize to larger molecules. 



Ef, + E]^ + Up — Ey — nu = 0, 



(10) 
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